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I . INTRODUCTION 


A quantitative  description  of  the  gas  flow  generated  by  a percussion 
primer  is  unavailable  at  the  present  time.*  In  this  report  we  describe 
a theoretical  model  and  calculational  results  of  the  flow  in  an  M-82 
percussion  primer  from  the  time  that  the  primer  propellant  is  ignited 
up  to  144  ys  after  the  gases  have  exited  the  vent  hole  and  enveloped 
an  inert  propellant  charge  in  the  chamber  of  the  gun. 

The  flow  down  the  flash  tube  from  the  primer  combustion  chamber  has 
been  modeled  as  one-dimensional,  time-dependent  flow  using  the  RECRIF 
code.  The  gas  motion  into  the  gun  chamber,  in  turn,  was  obtained  from 
the  revised  SAMS  code.^  Outside  of  the  primer  combustion  chamber  the 
flow  has  been  treated  as  that  of  a single-phase,  non-reactive  medium. 

In  Chapter  II,  we  describe  the  flow  model,  the  mathematical  aspects 
of  the  simulation  and  experimental  observations  used  in  the  validation 
of  the  results.  The  next  chapter  discusses  the  results  of  the  calcula- 
tions while  in  the  last  section  we  draw  general  conclusions  about  the 
results . 


II.  THE  FLOW  MODEL 


A.  The  Physical  Model 
1 . Overview 

_3 

The  M-82  percussion  primer  consist  of  22  grains  (1.43  x 10  kg) 
of .slow  burning  class  3 black  powder  occupying  a volume  of  1.528  x 10-6 
(nr)  of  the  forward  end  of  the  primer  combustion  chamber.  The  web  size 
of  the  propellant  was  taken  as  1.2  x 10  (m) . The  rest  of  the  chamber 

is  filled  with  an  inert  material.  The  combustion  chamber^is  connected 
by  a flash  tube  of  1.24  x 10  (m)  in  length  and  6.3  x 10  (m)  in  diame- 

ter through  an  opening,  called  the  vent  hole,  to  the  chamber  housing  the 
propellant  bags  (Figure  1) .+ 


* An  understanding  of  the  dynamical  processes  at  the  primer  charge 
interface  is  of  interest  to  charge  designers  in  their  effort  to 
achieve  smooth  pressure  transient  in  new  charge  configurations . 


Celmins , A.  K.  R.,  "RECRIF  Users'  Manual",  BRL  MR  2693,  October  1976. 

^Farr,  J.  L.  and  Traci,  R.  M. , "A  Users'  Manual  for  SAMS",  BRL  Contract 
Report  162,  June  1974. 

^See  Appendix  for  figures. 
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Upon  ignition, the  pressure  in  the  primer  chamber  increases.  A 
pressure  wave,  followed  by  the  combustion  products  of  the  propellant 
which  consist  of  burnt  gas,  unburnt  propellant  grains  and  solid  products 
of  combustion,  moves  down  the  tube,  finally  exiting  at  the  vent  hole.  To 
prevent  the  propellant  grains  from  blowing  out  of  the  combustion  chamber, 
a blowout  disc,  1.0  x 10”  (kg)  in  mass,  is  located  at  the  entrance  of  the 
flash  tub^  at  time  t = 0.  The  hold-back  pressure  of  the  disc  was  set  at 
2.07  x 10  (Pa)  ensuring  a pressure  buildup  in  the  combustion  chamber 
prior  to  the  disc  motioij.  The  bore  resistance  of  the  flash  tube  was 
assumed  to  be  1.67  x 10  (N)  and  the  friction  coefficient  assumed  to  be 

5.0  x 10  . 

The  vent  hole  conditions,  after  a short  transient,  correspond  to 
that  of  choked  flow  leading  to  an  underexpanded  jet  in  the  chamber  hous- 
ing the  propellant  bags. 

With  a stand-off,  that  is  a space  between  the  vent  hole  and  the 
propellant  bag,  as  in  the  present  study,  the  gas  upon  exiting  the  flash 
tube  expands.  The  flow  is  preceded  by  a weak  shock  which  impinges  upon 
and  is  reflected  from  the  base  of  the  propellant  bag.  The  gas  in  turn, 
due  to  the  geometrical  constraints  and  in  view  of  the  pressure  gradients 
present,  flows  along  and  around  the  propellant  bag,  eventually  complete- 
ly enveloping  it. 

-3 

In  this  study,  the  vent  hole_natthe  opening  was  6.35  x 10  (m)  in 

diameter,  the  stand-off  5.08  x 10  (m) , the  gun  tube  chamber  17.78  x 

10“2(m)  in  diameter  and  the  spacing  between  the  symmetrically  placed 
propellant  bag  and  the  gun  tube  wallwas  0.63  xl0"^(m). 

The  propelling  charge  of  the  round  consists  of  bags  of  propellant 
grains,  rather  closely  packed,  with  typically  a center  core  ignitor. 

The  breech  end  of  the  bag  is  terminated  by  a pad  of  black  powder  which 
is  ignited  by  the  hot  gases  of  the  primer,  thus  ultimately  leading  to 
the  ignition  of  the  whole  charge.  By  this  time  the  propellant  bag 
breaks  up,  the  propellant  grains  fill  the  whole  chamber,  resulting  in  a 
fluidized  mixture  of  propellant  and  combustion  gases.  With  the  increase 
in  pressure,  the  bore  resistance  is  overcome,  the  projectile  starts  to 
accelerate  down  the  tube,  enlarging  the  chamber  volume  available  for 
combustion  and  ultimately,  when  the  increase  in  volume  is  no  longer  com- 
pensated by  the  rate  of  gas  generation,  a decrease  in  chamber  pressure 
can  be  observed. 

An  accurate  tailoring  of  the  heat  and  mass  source  represented  by 
the  primer  is  then  essential  for  a smooth  ignition  of  the  propelling 
charge  of  the  gun. 
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2.  Flow  Through  the  Initiator 


Upon  initiation,  the  combustion  of  the  primer  grains  generates  a 
mass  of  high  pressure  gases*  These  gases  flow  down  a duct,  called  the 
flash  tube,  exiting  into  the  chamber  holding  the  propellant  bags.  For 
purposes  of  this  simulation,  it  was  assumed  that  the  combustion  chamber 
of  the  primer  is  connected  by  a constant  area  tube,  terminating  at  the 
vent  tube  to  the  stand-off  region.  Thus,  wave  reflection  due  to  changes 
in  geometry  is  neglected.  Only  product  gases  were  allowed  to  exit  the 
combustion  chamber  and  the  blow  out  disc,  after  it  reaches  the  vent  hole, 
is  automatically  removed  from  the  system,  i.e.,  no  longer  considered  as 
an  influence  on  the  flow  development. 

The  flow  in  the  flash  tube  is  assumed  to  be  one-dimensional,  lami- 
nar, compressible,  inviscid,  and  time-dependent;  and  chemical  reactions 
are  assumed  to  be  completed  by  the  time  they  exit  the  primer  combustion 
chamber.  The  details  of  the  chemical  kinetics  in  the  flow  are  not  known 
at  the  present  time,  but  by  using  a high  enthalpy  gas,  with  a temperature 
in  the  combustion  chamber  corresponding  to  that  of  the  flame  temperature 
of  the  propellant,  a good  approximation  of  the  flow  is  obtained.  The 
Noble-Abel  equation  of  state,  with  the  co-volume  taken  as  zero  was  used 
in  the  calculations.  In  Table  I,  the  pertinent  geometrical  and  propel- 
lant parameters  of  the  M-82  primer  are  listed. 


Table  I.  M-82  Primer  Characteristics 


Geometry 

Volume  occupied  by  the  Class  3 Black  Powder 
Flash  Tube  Diameter 
Flash  Tube  Exit  Diameter 
Flash  Tube  Length 
Combustion  Chamber  Volume 
Propellant  (Class  3 Black  Powder) 


9 . 310  x 10  (m  ) 


1.528  x 10  ^(m3) 
6.3  x 10~^(m) 
0.635  x 10"2(m) 


0.124 


(m) 

-6,  3 


Mass 

Burning  Rate  Exponent 
Burning  Rate  Coefficient 
Ratio  of  Specific  Heats 
Web  Size 
Explosive  Force 
Flame  Temperature 


1.43  x 10  ^(kg) 
7.540  x 10"1 
6.110  x 10"' 


1.20 


0.0012  (m) 
2.800  x 104  (m) 
2.030  x 103  (k) 
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To  account  for  the  presence  of  particles  in  the  flow  without  the 
benefit  of  a multiphase  analysis,  both  the  molecular  weight  and  the 
ratio  of  specific  heats  have  to  be  adjusted.  The  degree  of  modification 
needed  depends  on  the  mass  of  the  suspended  particles.  For  low  particle 
loading  satisfactory  flow  predictions  based  on  a one-phase  model  may  be 
obtained. 

Particle  laden  flows  have  been  studied  intensely  in  connection  with 
the  rocket  program  where  considerable  degradation  of  nozzle  performance 
was  observed. ^ In  these  studies,  particles  represented  up  to  30%  of  the 
total  mass  flow.  Only  rudimentary  studies  of  gas  particle  jets  exist. 
Two-phase  jets  have  been  treated  as  gas  jets  with  an  increased  density. 
For  a particle  jet  to  be  dynamically  equivalent  to  a gas  jet  the  gas  jet 
must  have  not  only  the  same  average  density,  but  must  also  exhibit  the 
same  mass  and  momentum  flux.  Injected  particles  generally  have  to  travel 
considerable  distances  before  they  reach  the  gas  velocity,  and,  as  Rudin- 
ger°  has  shown,  for  most  practical  cases,  this  is  an  unlikely  event. 

Neglecting  the  effect  of  turbulence,  assuming  uniform  distribution 
of  the  particles  over  the  cross  section  of  the  duct,  uniformity  in  size 
and  physical  properties,  precluding  evaporation,  condensation,  chemical 
reaction  or  breakup  of  the  particles  and  with  particle  size  of  the  order 
of  1 ym, the  analysis  is  straightforward. 

Within  the  context  of  a one-phase  flow  analysis,  one  can  account 
for  some  two-phase  flow  effects,  such  as  the  decrease  of  the  specific 
heat  ratio,  with  an  increase  in  mass  fraction  of  the  mixture. 

For  mass  fractions  as  low  as  0.03  and  a density  ratio  of 

£ ~ 0(10"3) , 
pp 

such  as  in  the  present  case,  the  particles  can  be  looked  upon  as  a spe- 
cies of  heavy  molecules,  with  a molecular  weight  that  is  several  orders 
of  magnitude  larger  than  the  gas.  In  effect,  one  assumes  then  that  the 
contribution  of  the  particles  to  the  pressure  is  negligible.  The  equa- 
tion of  state  for  the  mixture  then  follows  from  the  perfect  gas  law  but 
with  the  density  of  the  mixture.  As  long  as  the  particle  volume  can  be 
neglected,  a gas-particle  mixture  may  be  treated  as  a perfect  gas  with 
modified  thermal  properties. 


^Rudinger,  G.,  ’’Flow  of  Solid  Particles  in  Gases”,  AGARDograph  No.  222, 
AGARD,  Neuilly  Sur  Seine,  France,  1976. 
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Let  o - concentration  of  the  particles  and  a be  the  gas  concentra- 
tion.  The  mass  fraction  of  theparticles  is  given  by  4>  = + CT  > where 

p 

the  mixture  density  is  p = a + cr  • 

m p 

One  can  now  derive  a relationship  for  the  equation  giving  the  ratio 
of  specific  heats  for  the  gas  particle  mixture.  Following  Rudinger  y y 
the  internal  energy  of  the  mixture  can  be  expressed  as 


e = (l-d>)c  T + 4> ct 
m v Y v 

where  (1-<J>)  represents  the  mass  of  the  gas  and  x is  the  temperature  of 
the  particles  and  <J>  is  the  particle  mass  fraction.  The  enthalpy  can 
then  be  written  as 

h = e + E.  = (1-40  c T + d>c  r + 2. 

m m Pm  v Pm 

Recalling  the  definitions  of  specific  heat  and  substituting 


c 


v 

m 


(l-<JOcv  + 4>c 


( 1— 4> ) Cp 


+ 4>c 


consequently  Y , the  ratio  of  specific  heats  of  the  mixture  is  defined 
by  m 

Y = OHOCp  + 4>c 

m (1-4>D  7 4>c 

or  if  y = ratio  of  specific  heats  of  the  gas,  6 = — , where  c is  the 

P 

specific  heat  of  the  solid  particle  and  n the  ratio  of  the  mass  of  par- 
ticles to  the  mass  of  the  gas,  then 


4Rudinger , G.,  "Some  Effects  of  Finite  Particle  Volume  on  the  Dynamics 
of  Gas  Particle  Mixtures",  AIAA  Journal  3,  1217-1222,  (1965). 

^Rudinger,  G.,  "Relaxation  in  Gas-Particle  Flow"  in  Nonequilibrium  Flows, 
Part  I,  P.  P.  Wegener,  Ed.,  Marcel  Dekker , NY,  1969. 
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1 + r\6 

~ ^ 1 + y r\6  . 

(Note  that  6 % 0(1)  for  many  gas  particle  combinations  of  interest) . 

There  is  considerable  disagreement  in  the  literature  on  the  proper 
value  for  the  specific  heat  of  black  powder.  While  the  Encyclopedia  of 

Explosives^*  give  a value  of  0 . 2^^—j-^gJ  calculations  based  on  the  actual 

composition  of  black  powder,  assuming  75%  KNO^,  15%  charcoal  and  10-6 

sulphur  yields  a value  of  c = 0 . ^^gm-dteg^}'  This  value  is  valid  in  the 

temperature  range  from  300  to  700 [ K 3 and  represents  an  average  value. 

Using  6 = 2.080,  n = 0.03,  the  value  of  Ym  turns  out  to  be  1.22. 

At  n = 0.1  (10%  of  the  propellant  unburnt)  Ym  = 1-192.  Thus,  it  was 

felt  that  y =1.20  was  not  an  unreasonable  value  to  use  in  the  calcu- 
’m 

lations . 

For  a calculation  of  the  temperature  of  the  product  gases  the  molec- 
ular weight  has  to  be  known.  If  one  assumes  that  the  exhaust  gases 
contain  solid  particles  as  products  of  combustion  , which  appears  to  be 

the  case,  the  molecular  weight  is  6 1 * ^ [jTcmoXj  * 0n  t*ie  ot^er  ^ 

a solid  phase  doesn’t  precipitate,  the  molecular  weight  turns  out  to  be 

35.0  rA-1  Based  on  the  available  evidence,  we  felt  that  the  higher 
|kmolJ 

value  is  the  more  reasonable  value  to  use  in  this  simulation. 


3 . Flow  at  the  Base  of  the  Propellant  Bag 

The  underexpanded  jet  formed  at  the  vent  hole  spreads  laterally  and 
longitudinally,  eventually  contacting  the  base  of  the  propellant  bag,  at 
which  time  the  flow  structure  changes.  A shock  stands  off  the  surface, 
of  the  bag  and  a system  of  weak  shocks  may  be  observed  as  one  goes  radi- 
ally along  the  propellant  bag.  The  shock  enables  the  flow  to  equilibrate 
with  the  ambient  gas,  given  sufficient  room  for  lateral  expansion. 


6Federoff,  B.  T.  and  Sheffield,  0.  E.,  ’’Encyclopedia  of  Explosives  and 
Related  Items”,  Volume  2.  Picatinny  Arsenal,  Dover,  NJ , 1962, 

7 

Freedman,  E.,  Private  Communication. 
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In  reality,  due  to  the  hole  punched  in  the  propellant  bag  by  the 
primer  jet,  and  the  breakup  of  the  bag,  the  flow  is  considerably  more 
complicated,  but  neither  experimental  data  nor  a theoretical  study  of 
the  process  is  available  at  the  present  time.  For  the  case  at  hand,  it 
was  assumed  that  the  bag  retains  its  structural  integrity;  i.e.,  it  does 
not  move  or  burn  and  the  primer  flow  may  be  treated  analogously  to  a jet 
impinging  normally  onto  a flat  plate. 

The  jet  was  assumed  to  be  an  inviscid,  one-component  high- enthalpy 
but  nonreactive  gas.  The  confining  surfaces,  such  as  walls,  for  the  pur- 
pose of  this  investigation  were  assumed  to  be  inert  and  at  ambient  tem- 
perature. Also,  a no-mass  flux  condition  was  strictly  enforced  on  all 
surfaces . 

In  the  stand-off  space,  for  the  calculations  discussed  in  Chapter 
III,  neither  a Mach  disk  nor  the  stand-off  shock  could  be  detected.  These 
are  postulated  in  the  classical  theory  for  steady  processes  while  the 
case  at  hand  is  unsteady.  A quick  calculation  shows  that  if  we  assume 
steady  conditions  at  the  nozzle  corresponding  to  that  which  prevails  at 
t = 0,  both  the  stand-off  shock  and  the  Mach  disk  would  fall  outside 
the  physical  region  between  the  vent  hole  and  the  base  of  the  inert  bag. 
Thus,  it  appears  fair  to  suppose  that  the  flow  will  be  subjected  to  a 
number  of  wave  propagation  processes  and  at  no  time  will  an  equilibrium 
configuration  be  achieved. 

According  to  observations  of  the  impingement  of  non-uniform  axi- 
symmetric,  supersonic,  steady  jets  on  a perpendicular  flat  plate  in  the 
Mach  number  regime  1.4g<M<2.83,  a symmetrical  shock,  convex  toward  the 
nozzle,  spans  the  jet.  The  flow  beneath  the  impingement  shock  is  sub- 
sonic and  is  turned  outward  due  to  the  plate.  The  streamlines  are  de- 
flected outward  and  as  the  flow  moves  radially,  a reacceleration  with  a 
sonic  line  near  the  jet  edge  can  be  observed.  The  constant  pressure  up- 
per boundary  of  the  wall  jet  region  is  the  jet  edge  streamline  which  has 
undergone  a centered  expansion  upon  passing  through  the  shock.  Bubbles 
result  from  the  interaction  of  the  plate  shock  with  very  weak  shock  waves 
which  are  present  in  the  flow. 

The  significance  of  these  observations  is  that  there  are  circumstan- 
ces under  which  the  point  of  maximum  pressure  does  not  occur  on  the  stag- 
nation streamline  but  rather  off  the  center  line.  This  pressure  distri- 
bution has,  as  a consequence,  the  presence  of  a slowly  recirculating 
fluid  in  the  shock  layer  with  the  diameter  of  the  bubble  ranging  up  to 
0.8  times  the  jet  diameter.  These  flow  anomalies  may  have  some  effect 


g 

Kalghatgi,  G.  T.  , Hunt,  B.  L. , "The  Occurence  of  Stagnation  Bubbles  in 
Supersonic  Jet  Impingement  Flows",  Aeronautical  Quarterly,  27,  169-18S, 
1976. 
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upon  the  site  and  distribution  of  the  ignition  spots  on  the  propellant 
bag.  In  the  present  simulation  only  tenuous,  circumstantial  evidence 
for  these  mechanisms  was  observed.  However,  for  primers  with  larger 
stand-off  distances  calculated  by  higher  order  techniques  which  intro- 
duce less  smearing  than  present  in  this  work,  such  shock  phenomena  may 
possibly  be  observed. 

B.  Mathematical  Aspects  of  the  Simulation 

The  calculation  of  the  flow  from  the  combustion  chamber  through 
the  flash  tube  of  the  initiator-spindle  assembly  was  based  on  the  one- 
dimensional form  of  the  conservation  laws : 


(pA)t  + (puA)x  = 0 

(puA)t  + (u2pA)x  + Apx  + ApF  = 0 

[pA(e  + j u2)]t  + [puA(e  +|u“)  +2-  ]x~  qpA  = 0 

and  the  Noble-Abel  equation-of-state 

n = R T P 
p M 1-np 


with  n set  equal  to  zero.  These  equations  were  put  into  characteristic 
form  and  solved  numerically . ^ 


Space- averaged  values  of  the  variables  were  used  ^or  the  gas  in  the 
combustion  chamber.  The  propellant  obeyed  the  regression  law 


where  z(t)  represents  the  distance  that  the  propellant  has  burnt,  n,  3 
are  constants,  and  p is  the  pressure.  With  the  assumption  that  the 
kinetic  energy  of  the  gas  in  the  chamber  is  small  compcired  to  its  inter- 
nal energy,  the  average  pressure  can  be  approximated  by  the  stagnation 
pressure.  Thus,  using  an  energy  balance  of  the  chamber  with  the  burning 
rate  equation  and  the  continuity  equation,  all  the  flow  variables  can  be 
computed . 

The  two-dimensional,  time-dependent  part  of  the  calculation  in  the 
stand-off  region  was  based  on  the  Euler  equations.  The  finite  difference 
analog  of  the  conservation  laws  was  solved  by  a first-order,  two-step 
algorithm  which  has  been  described  elsewhere. ^ In  the  course  of  the 
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calculation,  the  shocks  are  captured  and  smeared  over  several  cell  widths, 
requiring  fine  gridding  for  the  resolution  of  flow  details  in  the  vicin- 
ity of  steep  gradients. 


A Noble-^bel  equation-of-state  was  used  with  the  co-volume  set  at 

1.17  x 10  ' [~]  which  is  valid  for  most  of  the  commonly  used  propellants. 

Points  along  the  solid  boundaries  were  calculated  by  the  reflection 
technique . 


C.  Experimental  Observations 

Only  a limited  amount  of  experimental  data  is  available  on  the 

primer  flow  dynamics. ^ Pressure  transducers  along  the  chamber 
shaft  of  the  initiator-spindle  assembly,  on  the  spindle,  as  well  as  on 
the  test  fixture  and  the  breech  end  of  the  propellant  bag,  yield  complete 
pressure  time  histories  of  the  product  gases  (see  Figure  1).  Further, 
using  high  resolution  photography,  the  expansion  of  the  gases  from  the 
vent  hole  into  the  test  fixture  without  a propellant  bag  have  been  ob- 
served (Figure  2) . From  these  records  the  velocity  of  the  gas  cloud  can 
be  deduced  and  compared  with  calculational  predictions.  Some  temperature 
measurements  have  also  been  attempted;  but,  at  present,  the  observations 

9 10 

are  felt  to  be  unreliable.  5 There  is  an  urgent  need  for  more  experi- 
mental data  on  all  facets  of  the  primer  flow  dynamics. 

9 

Stiefel  reports  that,  for  a gage  located  0.0508  [m]  from  the  axis 
of  the  tube  and  for  an  M-82  primer  of  the  M-198  howitzer  on  the  base  of 
the  propellant  bag  (P^,  Figure  1)  with  a stand-off  of  0.025  [m]  , the 
height  of  the  first  pressure  peak  varied  between  39  and  36  psi  (2.68  x 

10^ [Pa]  and  2.48  x 10^ [Pa]),  measured  at  1.2  and  2.8  [ms]  respectively 
for  a straight-through  duct.  With  a squared  cone  exit  the  pressure 
peaks,  though  of  the  same  magnitude,  were  observed  earlier,  i.e.,  at 
0.6  [ms].  However,  in  changing  from  a straight-through  hole  to  a squared 
cone  exit,  the  pressures  exhibited  more  scatter;  that  is,  they  ranged 

between  33  and  44  psi  (2.27  x 10^[Pa]  and  3.03  x 10^[Pa])  in  the  experi- 
mental runs. 


Stiefel,  L.,  "Experimental  Characterization  of  the  Output  of  Artillery 
Primers",  (to  be  published). 

^White,  K.  J.,  Hartman,  R.  A.,  May,  I.  W. , Kelso,  J.  R.,  "Experimental 
Investigation  of  Ignition  Train  System  for  Bagged  Charges",  14th  JANNAF 
Combustion  Meeting,  (in  press),  1977. 

^White,  J.  J.,  Price,  C.  F.,  Holmes,  H.  E.,  May,  I.  W. , "Black  Powder 
and  Clean  Burning  Ignitor  Train  Studies",  13th  JANNAF  Combustion  Meet- 
ing, Vol . I,  pp.  405-423,  CPIA,  Laurel,  MD,  1976. 
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The  experiments  further  showed  that  the  velocity  of  the  pressure 

m 9 

wave  is  around  365  [— ] in  the  axial  direction,  while  laterally,  Stiefel 

estimated  a velocity  of  40  [— ] . White^,  based  on  time  lapse  photogra- 
phy of  a primer  jet  exhausting  into  a chamber,  found  that  the  center 
core  of  the  jet,  at  a distance  of  0.1  [m]  from  the  vent  hole, propagates 

at  300  [-]. 

Ls J 


III.  RESULTS 

A.  Flash  Tube  Calculations 

The  calculat ional  results  can  be  divided  into  two  parts:  the 
flow  in  the  flash  tube  between  the  primer  combustion  chamber  and  the 
vent  hole  and  the  flow  in  stand-off  space  between  the  vent  hole  and  the 
propellant  bag. 

The  one-dimensional  calculations  of  the  flow  in  the  primer  combus- 
tion chamber/f lash  tube  were  carried  out  using  the  RFCRIF  code.  First, 
the  burning  rate  exponent  and  coefficient  were  obtained  from  experimen- 
tal data.  These  parameters  and  y,  the  ratio  of  specific  heats,  were  var- 
ied until  the  pressure-time  behavior  at  the  station  0.084 [m]  from  the 
combustion  chamber  of  the  primer  could  be  matched  at  the  high  end  of  the 

-7 

pressure-time  curve  (see  figure  3).  A value  of  8 of  6.10  x 10  with  n, 
the  burning  rate  exponent  equal  to  0.754  and  a y = 1.20  gave  acceptable 
results.  The  calculation  of  the  flow  over  the  whole  length  of  the  flash 
tube  followed,  with  our  primary  interest  centering  on  the  vent  hole 
conditions . 

-4 

At  t = 5.00  x 10  [s]  after  the  ignition  of  the  black  powder,  the 

propellant  gases  reach  the  vent  hole  and,  shortly  afterwards,  the  flow 
is  choked.  In  Figure  4 we  have  a plot  of  the  mass  of  combustion  products 
as  a function  of  time  in  seconds.  The  uppermost  curve  gives  the  total 
mass  burnt,  while  the  lowest  curve  represents  the  propellant  gas  mass 
which  has  exited  the  flash  tube.  The  difference  in  the  ordinate  for  each 
time  between  the  two  lower  curves  gives  the  mass  which  is  in  the  flash 
tube,  while  the  difference  between  the  two  topmost  curves  gives  the 
amount  of  mass  still  residing  in  the  combustion  chamber.  Since  we  know 
the  total  propellant  mass  initially  available,  one  can  estimate  the 

.3 

unburnt  mass  at  each  time.  For  example,  if  at  t = 1.514  x 10  [s]  we 

assume  that  the  unburnt  particles  are  all  evenly  distributed  within  the 
flash  tube  and  chamber,  one  arrives  at  a figure  of  3%  by  mass  of  solid 
particles  in  the  flow. 

The  velocity-time  history  of  the  flow  in  the  flash  tube  is  illustra- 
ted in  Figure  5.  The  lower  curve  gives  the  velocity  at  the  entrance  to 
the  flash  tube  from  the  primer  combustion  chamber  while  the  other  curve 
gives  the  velocity  at  the  vent  hole.  This  consists  of  two  parts: 
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the  rising  portion  gives  the  velocity  of  the  blowout  disc  (i.e.,  the 
-3 

1.0  x 10  [kg]  Plug  which  was  used  to  prevent  the  unburnt  particles  from 
ejecting  from  the  primer  duct),  while  the  decreasing  portion  gives  the 
velocity  at  the  muzzle.  It  should  be  noted  that  initially  the  flow  is 
supersonic  at  the  muzzle,  but  it  drops  continuously  until  sonic  condi- 
tions are  obtained.  M = 1 occurs  at  the  break  in  the  slope  of  the  curve. 

The  next  two  curves  give  the  vent  hole  temperature,  in  degrees  [K] 
and  the  vent  hole  pressure  in  [Pa]  as  a function  of  time  [s].  Zero  time 
refers  to  the  inception  of  burning  in  the  combustion  chamber  and  the 
steep  jump  occurs  as  the  flow  reaches  the  vent  hole.  The  break  in  the 
slope  of  the  pressure  curve  again  indicates  transition  to  the  choked 
condition . 

From  the  data  of  the  pressure,  density  and  velocity  at  the  vent  hole, 
one  obtains  the  following  analytical  expressions  for  these  quantities  as 
a function  of  time: 

P = 0.128  x 1014  x t2  - 0.302  x 1011  x t + 2.140  x 10? 

p = 0.186  x 108  x t2  - 0.791  x 105  x t + 9.141  x 10 

u = 0.802  x 108  x t“  - 1.451  x 105  x t + 4.877  x 102 

The  next  four  figures  (8  through  11)  give  the  velocity,  temperature, 

pressure  and  density  in  the  chamber  and  the  flash  tube  at  four  different 

times,  i.e.,  at  t = 3.733  x 10"4,  7.466  x 10_4,  1.120  x 10'3  and  1.493  x 

“ 3 

10  [s].  The  vertical  line  at  the  left  is  the  separation  point  between 
the  combustion  chamber  and  the  flash  tube.  The  velocity  appears  as  a 
linearly  increasing  function  of  distance  while  the  temperature,  except 

-4 

for  the  initial  time,  t = 3.733  x 10  [s],  is  a slowly  decreasing  func- 

tion. It  stays  below  the  flame  temperature  since, as  the  gases  expand, 
cooling  takes  place.  This  type  behavior  is  as  expected.  The  initial 

-2  . 

curve  terminates  at  5.9  x 10  ^[m],  i.e.,  at  the  location  of  the  blowout 
disc  within  the  flash  tube.  At  succeeding  times  in  these  figures,  the 
disc  has  already  exited  the  tube.  The  lower  three  curves  terminate  at 
the  vent  hole  exit. 

The  gas  velcjpity  at  the  venthole,  as  shown  in  Figure  8,  is  at  a 
velocity  of  420  [— ] , but  as  time  progresses  and  the  pressure  drops,  it 
too  drops.  It  is  interesting  to  note  the  almost  constant  pressures  and 
densities,  at  each  time,  all  along  the  flash  tube. 
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The  conditions  at  the  vent  hole,  which  are  the  input  conditions  to 
the  second  phase  of  the  calculations,  that  is  the  two-dimensional  time- 
dependent  flow  in  the  stand-off  zone,  were  followed  for  1.514  [ms]. 

The  pressure,  velocity  and  density  were  then  curve-fitted,  using  a Lagran- 
gian  interpolation  technique,  to  give  an  analytical  expression  for  the 
time  behavior  of  these  functions . The  formulas  were  checked  at  various 
points  along  the  time  axis  against  the  calculated  values  to  insure  the 
range  of  validity  of  the  interpolation  technique. 

B.  Flow  in  the  Stand-off  Region 


Once  the  gases  exit  the  vent  hole,  their  behavior  changes  dras- 
tically. A one-dimensional  analysis  is  no  longer  adequate  as  it  is 
not  able  to  account  for  the  effect  of  the  lateral  motion  in  the  space 
between  the  base  pad  and  the  vent  hole.  The  computational  results  ob- 
tained by  exercising  the  revised  SAMS  code  will  now  be  summarized. 

With  the  time-dependent  vent  hole  conditions  as  input,  the  initial 
runs  calculated  a free  jet  formed  by  the  gas  issuing  from  the  flash  tube. 
The  results  were  then  compared  with  the  photographs  shown  in  Figure  2(a) 
through  2(d). 

Figures  12  and  13  show  pressure  contour  plots  of  the  jet  with  no 
charge  in  place  at  23.14  and  66.19  [ys] , respectively.  The  vent  hole  is 
located  in  the  upper  left-hand  corner  and  the  axis  of  the  jet  is  at  the 
left-hand  border  of  the  plot.  At  t = 23.14  [ys],  for  example,  the  pres- 

P . 

sure  ratio  across  the  plane  of  the  vent  hole  is  p^=  1 . 5,  resulting  in  an 

underexpanded  jet.  The  closely  spaced  contour  lines  indicate  the  pres- 
ence of  a shock.  As  time  progresses  the  shock  weakens  and  also  diffuses. 
This  can  be  seen  very  clearly  in  the  pressure-distance  plot  (along  the 
jet  axis)  for  the  three  different  lines  (Figure  14)  where  the  pressure 
hump  has  spread  and  its  amplitude  lowered  when  the  two  curves  are 
compared.  The  density  plots,  Figures  15  and  16,  exhibit  a beha\ior  anal- 
ogous to  that  of  the  pressure. 

_2 

Along  the  spindle,  in  the  vicinity  of  2.0  x 10  [m]  from  the  center- 

line  of  the  jet,  one  sees  a region  of  lower  pressure,  which  is  followed 
by  a steep  pressure  rise  as  one  approaches  the  blast  wave. 

The  radial  pressure  distribution  across  the  jet  at  three  different 
times  is  shown  in  Figure  17,  while  near  the  tip  of  the  jet,  the  pressure 
hump  due  to  the  blast  wave  is  clearly  discernible  (Figure  18). 

The  velocity  in  the  axial  direction  along  the  jet  for  t = 24.14, 
45.73,  and  66.19  [ys]  is  shown  in  Figure  19.  Contrary  to  the  contour 
plots,  one  should  note  that  in  this  and  the  succeeding  figures,  the 
origin  of  the  x-axis  is  at  the  vent  hole  and  the  coordinate  increases 
in  the  downstream  direction. 
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As  the  gas  leaves  the  flash  tube,  the  pressure  ratio  at  the  vent  hole  at 
the  latter  time  is  1.06,  so  the  conditions  for  an  underexpanded  jet  are 
still  clearly  met.  The  gases  expand  and  accelerate.  A blast  wave  is 

seen  to  move  with  a peak  velocity,  at  66.19  [us],  of  around  1500  [ j]  . 

The  main  body  of  the  gas  moves  slower;  near  the  center  of  the  plume  it 

is  around  600  [— ] . 

Ls 

From  the  photographs  of  Figure  2 one  sees  that  the  center  of  the 

plume  propagates  at  around  300  [— ] . This  points  out  the  main  shortcom- 
ing of  the  present  simulation.  With  a code  which  is  restricted  to  a 
single  y for  both  the  high  enthalpy  reaction  products  and  the  ambient, 
some  of  the  calculated  flow  parameters  will  be  inaccurate.  We  set  the 
ambient  density  at  the  value  corresponding  to  air  at  300  [K]  and,  conse- 
quently, the  ambient  pressure,  due  to  the  y used,  becomes  0.485  x 10  [Pa] 
instead  of  the  1.013  x 10b[Pa].  Elementary  shock  wave  theory  then  shows 
that  the  rate  of  shock  propagation  will  be  too  high  by  a factor  of  a little 
over  2.  On  the  other  hand,  the  peak  pressures,  as  shown  below,  can  be 
calculated  quite  accurately. 

With  the  propellant  bag  in  place,  modeled  as  a reflecting  wall  at 
0.025[m]  from  the  vent  hole,  the  jet  cannot  proceed  far  downstream  be- 
fore it  is  turned  radially.  The  pressure  contour  plots,  Figures  20  and 
21,  give  an  overall  view  of  the  development  of  the  flow.  Again,  at  the 
edge  of  the  plume  a strong  blast  is  in  evidence  which  weakens  as  it  moves 
away  from  the  jet  axis.  A region  of  low  pressure  behind  the  blast  is  in 
evidence  in  Figure  20b  and  even  more  clearly  in  Figure  21.  At  t = 144.40 
[ys],  near  the  end  of  the  simulation,  the  blast  wave  has  weakened  to  the 
point  where  the  pressure  ratio  at  the  head  of  the  wave,  with  respect  to 
the  conditions  ahead  of  it,  is  less  than  3.  The  density  contour  plots, 
Figures  22  and  23,  mimic  those  of  the  pressure. 

Considerable  computational  information  is  avaiable  for  the  pressure 
distribution  within  the  flow  field  and  comparison  with  experimental  data 
is  possible.  Along  the  jet  axis  (Figures  24  through  26,  from  the  vent 
hole  to  the  propellant  bag,  the  pressure  is  decreasing  from  an  initial 

value  of  1.6  x 107[Pa]  at  t = 74.73[ys],  for  example,  at  the  vent  hole, 

to  1.05  x 105[Pa]  at  the  stagnation  point.  Over  the  time  range  of  this 
calculation,  as  the  vent  hole  pressure  changed  only  marginally  from  2.14 

x 107[Pa]  to  1.5  x 10' [Pa]  at  t = 144.40[ps],  quasi-steady  conditions 
were  observed  at  the  propellant  bag.  After  the  initial  transient,  the 
pressure  at  the  stagnation  point  changed  by  less  than  11%  so  that  condi- 
tions there  corresponded  to  a quasi -steady  state.  This  is  no  longer 
true  as  one  proceeds  radially  along  the  base  of  the  propellant  bag. 
In  Figure  27  one  can  clearly  see  the  blast  wave  sweeping  the 
surface.  The  flow  is  expanding  behind  the  wave  with  the  pressure 
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continuously  dropping.  At  t = 46.44[ys]  in  Figure  27,  the  lowest  pres-^ 
sure  behind  the  blast  is  5.0  x 10  [Pa]  at  r = 0.020  [m]  from  the  jet  axis 
It  drops  continuously  as  can  be  seen  in  the  next  two  figures.  In  Figure 
28,  the  first  indication  of  a pressure  plateau  is  in  evidence  which  be- 
comes very  pronounced  at  t = 144.40  [us]  in  Figure  29.  By  this  time  the 

pressure  behind  the  blast  had  dropped  to  1.0  x 10^ [Pa]. 

The  flow  exhibits  analogous  characteristics  if  one  traverses  it  in 

-3 

the  radial  direction  at  a distance  of  4.00  x 10  " [m]  from  the  propellant 
bag  base  surface.  Again,  as  shown  in  Figure  30,  after  an  initial  pres- 
sure ramp  at  early  times,  clearly  a blast  wave  develops  which  propagates 
rapidly  outward.  Later,  for  example  at  t = 102.63  [us]  in  Figure  31,  sub 
atmospheric  conditions  prevail  behind  the  blast  followed  by  a pressure 
rise  back  to  the  vent  hole  conditions  as  one  mo\res  toward  the  jet  axis. 

In  these  diagrams  the  pressure  ramp  is  a consequence  of  the  diffusion 
of  the  wave. 

The  next  series  of  diagrams,  Figures  33  through  35,  show  the  velo- 
city distribution  along  the  jet  axis.  Again,  after  an  initial  transient, 

the  maximum  velocity  stays  at  1200  [—]  . 

In  Figure  36  a comparison  is  shown  between  the  calculated  and  the 
observed  temporal  distribution  of  the  pressure  at  a location  of  0.05  [m] 
radially  from  the  tube  axis  at  the  base  of  the  propellant  bag  (Figure  1 

gage  P.) . Stiefel’s  data9  shows  a pressure  peak  of  2.48  x 105[Pa]  at 
0.6  [ms]  with  a squared  cone  exit.  For  a straight-through  hole,  the  peak 
pressure  was  observed  at  a later  time,  somewhere  between  1.2  and  1.6  [ms] 
On  the  other  hand,  the  peak  pressure  did  not  seem  to  be  influenced  appre- 
ciably by  the  geometry  of  the  vent  hole.  We  see  that  the  calculated  pres 
sure  follows  closely  the  observed  behavior.  The  only  difference  is  that 
after  100  [ys] , the  calculated  pressure  drops  off  more  steeply  than  the 
measured  value.  The  difference  between  the  calculated  and  the  experi- 
mental values  at  these  late  times  may  be  due  to  the  frequency  response 
of  the  gage. 

The  temperature  along  the  base  of  the  propellant  bag  measured  from 
the  centerline  of  the  bag  reflects  the  passing  of  the  blast  wave  and, 
later,  the  propellant  gases.  Starting  at  the  stagnation  point  and  mov- 
ing radially,  one  notices  the  first  local  temperature  maximum  at  around 

7 

5.0  x 10_o[m]  which  remains  in  the  same  location  from  t = 16.78  [us]  to 
the  termination  of  the  computer  run.  This  evidence  points  to  the  pre- 
sence of  a shock,  made  plausible  by  the  arguments  advanced  in  Chapter  II, 
Section  A-3,  of  this  report. 

The  second  local  temperature  maximum  is  observed  at  each  time  frame 
at  the  location  of  the  blast  wave  where  the  temperature  is  near  1100[Kl. 
The  temperature  data  for  several  times  is  summarized  below  in  Table  II. 
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Table  II.  Temperature  Along  the  Base  Pad 


Time  [ys] 

60 . 66 

102.63 

144.40 

Location  [m] 

T[k] 

T[K] 

t[k] 

0.001 

629.58 

617.18 

604.47 

0.005 

736.10 

721.52 

706.61 

0.010 

739.47 

724.86 

709.89 

0.015 

720.96 

706.75 

692.16 

0.020 

696.59 

682.88 

668.80 

0.025 

674.28 

656.51 

643.00 

0.030 

1045.58 

629.41 

616.39 

0.035 

665.05 

639.27 

590.40 

0.040 

477.18 

1070.07 

755.94 

0.045 

301.57 

814.12 

1015.38 

0.050 

480.78 

850.85 

0.055 

398.85 

544.92 

0.060 

367.85 

398.97 

0.065 

300.62 

357.66 

0.070 

322.29 

0.075 

340.20 

0.080 

300.61 

0.085 

300.00 

These  temperatures  should  not  be  considered  as  absolute  values  but 
rather  as  indicative  of  the  thermodynamic  state  at  that  particular  time 
and  place.  As  pointed  out  earlier,  the  temperature  is  a derived  quan- 
tity, obtained  from  the  equation-of-state . Depending  on  the  value  of  the 
molecular  weight  of  the  product  gases,  the  temperature  can  be  made  to 
vary  by  a factor  of  two.  The  controversial  point  is  the  question  whe- 
ther solids  are  precipitated  as  products  of  the  reaction,  and  if  so,  in 
what  quantities.  A judgment  of  the  accuracy  of  the  temperature  predic- 
tion has  to  await  the  resolution  of  this  kinetics  question. 
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It  should  be  kept  in  mind  that  all  of  the  foregoing  results  were 
obtained  under  greatly  simplifying  assumptions;  the  most  restrictive  of 
these  was  the  fact  that  only  one  species,  i.e.,  only  one  value  of y was 
used  for  both  the  reactants  and  the  ambient.  They  used  was  that  of 
the  high  enthalpy  gas.  Thus,  the  ambient,  though  at  the  correct  tem- 
perature, is  at  a pressure  of  5.22  x 10^ [Pa]  which  is  lower  than  atmos- 
pheric, leading  to  an  artificially  enhanced  wave  propagation  through 
the  medium. 

While  overall  determination  of  the  flow  dynamics  is  not  adversely 
affected  by  the  stipulation  of  laminar  flow  and  negligible  transport 
effects,  local  behavior  can  be  greatly  distorted  by  these  assumptions. 
The  inclusion  of  these  effects  would  require  the  solution  of  the  Navier- 
Stokes  equations  with  kinetics  and  a turbulence  model.  At  the  present 
time  the  interaction  between  turbulence  and  chemical  reactions  is  in- 
completely understood. 

12  13 

Turbulence  models  for  compressible  flows  * have  been  developed 
and  applied  with  some  degree  of  success  to  practical  problems.  However 
these  theories,  as  yet,  do  not  take  cognizance  of  the  presence  of  par- 
ticulate matter  in  the  flow. 

The  overriding  interest  in  primer  flow  calculations  is  the  possi- 
bility of  the  determination  of  the  heat  transfer  to  the  propellant  bed. 
Known  correlations  relating  Nusselt  to  the  Reynolds  and  Prandtl  numbers 
of  course  require  values  of  the  transport  coefficients  which  were  not 

14  15 

considered  in  this  study  * . However,  from  a study  of  the  flow  re- 

gime, one  can  make  reasonable  estimates  of  these  and  thus  of  the  heat 
transfer.  Ultimately,  of  course,  the  full  three-dimensional,  two-phase, 
unsteady,  viscid,  reactive  flow  equations  will  have  to  be  solved  for  a 
complete  resolution  of  the  problem. 
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IV.  CONCLUSIONS 


A theoretical  program  has  been  conducted  to  investigate  the  flow 
dynamics  of  an  M-82  primer  from  the  time  that  the  primer  propellant  is 
ignited  up  to  144  [ys]  after  the  gases  have  exited  the  vent  hole  and 
enveloped  an  inert  propellant  charge  in  the  chamber  of  the  gun.  The 
results  of  this  study  indicate: 

a.  The  gross  features  of  the  primer  flow  dynamics  can  be  modeled 
adequately  with  off-the-shelf  computer  codes; 

b.  These  results  permit  an  overall  assessment  of  primer  functioning; 

c.  Comparison  with  experimental  findings  is  only  possible  for  some 
of  the  parameters  of  the  problem; 

d.  Using  available  codes,  parametric  studies  can  be  conducted  which 
should  point  to  modes  of  product  improvement. 
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STAND  OFF  REGION 
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Figure  1.  Schematic  of  the  initiator  spindle  assembly  for  the  155mm  howitzer.  P and  P 
indicate  the  location  of  the  pressure  transducers  (Reference  g). 


Figure  2.  The  development  of  the  flow  at  the  vent  hole.  (a}  t = 0,  52.9  [us] ; 

(b)  t = 106.  159  [us]  ; (c)  t = 212,  265  [us];  (d)  t = 317,  370  [us]. 
Right-hand  side  of  photo  is  viewed  from  90°.  Each  small  division 
represents  1.0  x 10"  [n] . Photographs  courtesy  of  Dr.  K.  .T.  Wliite, 
Propulsion  Division,  RRL. 


52 


d 


Figure  2 
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Figure  3.  Comparison  of  calculated  and  experimentally  observed  pressure 
time  history  for  black  powder  at  a station  0.084  [m]  from  the 
primer  combustion  chamber.  Experimental  data  courtesy  of 
Dr.  K.  J.  White,  Propulsion  Division,  BRL. 
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Figure  4.  Mass  of  combustion  products  in  [kg],  as  a function  of  time  in  [s]  in  the  flash  tube. 
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Figure  5.  Gas  velocity  in  [— ] as  a function  of  time  in  [s]. 
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Figure  6.  Vent  hole  temperature  in  degrees  [K]  as  a function  of  time  expressed  in  [s]. 
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Figure  7.  Vent  hole  pressure  in  units  of  [Pa]  as  a function  of  time  in  [sj. 
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Figure  9.  Temperature  in  the  chamber  and  the  flash  tube  as  a function  of  location.  The  vertical 
line  denotes  the  entrance  to  the  flash  tube.  Units  are  [K]  and  [m]  respectively.  _ 
From  top  to  bottom  the  curves  refer  to  t = 3.733  x 10“^,  7.466  x 10”  , 1.120  x 10~ 
and  1.493  x 10~^[s]  respectively. 
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Figure  10.  Pressure  in  the  chamber  and  the  flash  tube  as  a function  of  location.  The  vertical 
line  denotes  the  entrance  to  the  flash  tube.  Units  are  [Pa]  and  [m]  respectively. 
From  top  to  bottom  the  curves  refer  to  t = 3.733  x 10“^,  7.466  x 10”^,  1.120  x 10“^ 
and  1.493  x 10“^[s]  respectively. 
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Figure  11.  Density  in  the  primer  combustion  chamber  and  the  flash  tube  as  a function  of  location. 

The  vertical  line  denotes  the  entrance  to  the  flash  tube.  Units  are  [kg/m3]  and  [m] 
respectively.  From  top  to  bottom  the  curves  refer  to  3.733  x 10-4,  7.466  x 10"4, 

1.120  x 10-3  and  1.493  x 10“ 3 [s]  respectively. 
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Figure  13.  Free  jet  pressure  contour  plot  at  t = 66.19  [ys] . 
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Figure  14.  Pressure  vs  distance  along  the  jet  axis. 
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Figure  16.  Free  jet  density  contour  plot  at  t = 66.19  [ys] . 
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Figure  17.  Pressure  radially  across  the  free  jet  at  1.9  x 10  [m]  from  the  plane  of  the  vent  hole 
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Figure  18.  Pressure  radially  across  the  free  jet  at  4.9  x 10  [m]  from  the  plane  of  the  vent  hole. 
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Figure  19.  Velocity  along  the  free  jet  axis  at  t = 24.14,  45.73  and  66.19  [ys]  respectively. 
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Figure  24.  Pressure  along  the  jet  axis  in  the  stand-off  region  at  t = 16.78 
31.97  and  46.44  [ys] . 
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Figure  25.  Pressure  along  the  jet  axis  in  the  stand-off  region  at  t = 60.66 
74.73  and  88.70  [us]. 
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Figure  27.  Pressure  radially  from  the  jet  axis  along  the  propellant  bag  base 
at  t = 16.78,  31.97  and  46.44  [ys] . 
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Figure  28.  Pressure  radially  from  the  jet  axis  along  the  propellant  bag  base 
at  t = 102.63,  116.54  and  130.46  [ys] . 
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Figure  29.  Pressure  radially  from  the  jet  axis  along  the  propellant  bag  base 
at  t = 144.40  [ys] . 
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Figure  30.  Pressure  radially  from  the  jet  axis  at  4.00  x 10  [m]  from  the  propellant 
bag  base  at  t = 16.78,  31.97  and  46.44  [ys] . 
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HORIZONTAL  VELOCITV 
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Figure  33.  Velocity  along  the  jet  axis  in  the  stand-off  region  at  t = 16.78, 
31.97  and  46.44  [ys] . 
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Figure  34.  Velocity  along  the  jet  axis  in  the  stand-off  region  at  t = 102.63, 
116.54  and  130.46  [ps]. 
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Figure  35.  Velocity  along  the  jet  axis  in  the  stand-off  region  at  t = 144.40  [vs]. 


PRESSURE  (Pa) 


106 


• CALCULATED 
— STIEFEL'S  DATA 


104L_i_J i I i . j 1 i 1 i I 

0 20  40  60  80  100  120  140 

TIME  (/is) 

Figure  36.  Comparison  between  calculated  and  measured  pressures  at  a 

location  0.0508  [m]  from  the  jet  axis  at  the  propellant  bag 
base  as  a function  of  time.  Stiefel's  data  is  from  Reference 

9. 
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LIST  OF  SYMBOLS 


A(x) 

base  area 

[m2] 

c 

specific  heat  of  solid  particles 

[J/K  kmol] 

c 

p 

specific  heat  at  constant  pressure 

[J/K  kmol] 

cv 

specific  heat  at  constant  volume 

[J/K  kmol] 

F 

force  per  unit  mass  due  to  friction 

[m/s2] 

P 

pressure 

[Pa] 

q 

heat  generation  rate/mass 

[J/kg*s] 

t 

time 

[s] 

u 

axial  gas  velocity 

[m/s] 

R 

universal  gas  constant 

[J/K  kmol] 

M 

molar  mass 

[kg/mol] 

n 

burning  rate  exponent 

z 

surface  regression  distance 

[a] 

8 

burning  rate  coefficient 

6 

c 

c 

P 

Y 

ratio  of  specific  heats  of  the  gas 

Ym 

ratio  of  specific  heats  of  the  mixture 

0 

density 

[kg/m3] 
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